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ABSTRACT 

We propose a physical model for nonlinear stochastic biasing of one-point statistics 
resulting from the formation epoch distribution of dark halos. In contrast to previous 
works on the basis of extensive numerical simulations, our model provides for the first 
time an analytic expression for the joint probability function. Specifically we derive 
the joint probability function of halo and mass density contrasts from the extended 
Press-Schechter theory. Since this function is derived in the framework of the standard 
gravitational instability theory assuming the random-Gaussianity of the primordial den- 
sity field alone, we expect that the basic features of the nonlinear and stochastic biasing 
predicted from our model are fairly generic. As representative examples, we compute 
the various biasing parameters in cold dark matter models as a function of a redshift 
and a smoothing length. Our major findings are (1) the biasing of the variance evolves 
strongly as redshift while its scale-dependence is generally weak and a simple linear bias- 
ing model provides a reasonable approximation roughly at R ^ 2(l-|-z)/i~^Mpc, and (2) 
the stochasticity exhibits moderate scale-dependence especially on R ^ 20/i~^Mpc, but 
is almost independent of z. Comparison with the previous numerical simulations shows 
good agreement with the above behavior, indicating that the nonlinear and stochastic 
nature of the halo biasing is essentially understood by taking account of the distribution 
of the halo mass and the formation epoch. 

Subject headings: cosmology: theory - galaxiesxlustering - galaxies:clusters:general - 
galaxiesiformation - galaxies: halos - dark matter - large-scale structure of universe 



1. Introduction 

The universe after the last scattering can be probed only through observing the distribution of 
luminous objects, either directly or indirectly via the weak lensing effect. This is why several wide- 
field and/or deep surveys of galaxies, clusters and quasars are planned and operating at various 



- 2 - 



wavelengths. The purpose of such cosmological surveys is two-fold; to understand the nature of 
the astronomical objects themselves, and to extract the cosmological information. From the latter 
point of view, which we will pursue throughout this paper, the objects serve merely as tracers of 
dark matter in the universe. Since such luminous objects should have been formed as a consequence 
of complicated astrophysical processes in addition to the purely gravitational interaction, it is quite 
unlikely that they faithfully trace either the spatial distribution of dark matter or its redshift 
evolution. Rather it is natural to assume that they sample the dark matter distribution in a biased 
manner. 

To describe the biasing more specifically, define the density contrasts of galaxies and mass at 
a position x and redshift z smoothed over the scale R as 

(5gal(x, Z\R) = 1, (^mass(x, z\R) = = 1. (1) 

'^gal Pmass 

Here and in what follows we use the words, "mass" and "dark matter", interchangeably, and 
"galaxies" to imply luminous objects (galaxies, clusters, quasars, etc.) in a collective sense. 

The simplest, albeit most unlikely, possibility is that they are proportional to each other: 

5gal (x, Z I i?) = 6 (^mass {^,z\R). (2) 

While the proportional coefficient was assumed to be constant when the concept of the biasing was 
first introduced in the cosmology community (Kaiser 1984; Davis et al. 1985; Bardeen et al.l986), 
it has been subsequently recognized that b should depend on z and R (e.g.. Fry 1996; Mo & While 
1996). As a matter of fact, it is more realistic to abandon the simple linear biasing ansatz (2) 
completely, and formulate the biasing in terms of the conditional probability function P(5gai|'^mass) 
of (5gai at a given dmass- Then equation (2) is rephrased as 

/oo 
<^gal -P(^gal|<^mass) dSggX- (3) 

Obviously the relation between (5gai and (5mass is neither linear nor deterministic. This general 
concept - nonlinear and stochastic biasing - was introduced and developed in a seminal paper 
by Dekel & Lahav (1999), which inspired numerous recent activities in this field (e.g.. Pen 1998; 
Taruya, Koyama & Soda 1999; Blanton et al. 1999; Matsubara 1999; Taruya 2000; Somerville et 
al. 2000; Taruya et al. 2000). 

Then the crucial question is the physical interpretation of P(^gail^mass)- There are (at least) 
two different interpretations for its physical origin. The first is based on the fact that Sg^i is 
meaningless unless one specifies many hidden variables characterizing the galaxies in the catalogue, 
for instance, their luminosity, mass of dark matter and gas, temperature, physical size, formation 
epoch and merging history, among others. This list is already long enough to convince one for the 
inevitably broad distribution of P(<5gai|5mass)- In this spirit, Blanton et al. (1999) proposed that 
the gas temperature of a local patch is the important variable to control -P((5gai|(5mass) on the basis 
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of cosmological hydrodynamical simulations. The other adopts the view that our universe is fully 
specified by the primordial density field of dark matter. According to this interpretation, all the 
properties of any galaxy should be in principle computable given an initial distribution of dark 
matter field in the entire universe, at least in a statistical sense. Actually this is exactly what the 
cosmological hydrodynamical simulations attempt to do. The gas temperature of a local patch, for 
instance, should be determined by a non-local attribute of the dark matter fluctuations. Clearly the 
above interpretations are not conflicting, but rather stress the two different aspects of the physics 
of galaxy formation which is poorly understood at best. 

In this paper, we present an analytical model for nonlinear stochastic biasing by combining the 
above two interpretations in a sense. Specifically we derive the joint probability function of 5haio 
and 5mass7 ^((^haio, ^mass)) from the distributions of the formation epoch and mass of halos on the 
basis of the extended Press-Schechter theory. In contrast to previous work which were based on the 
results of numerical simulations (Kravtsov k. Klypin 1999; Blanton et al. 1999; Somerville et al. 
2000), our model provides, for the first time, an analytic expression for -P((5haio5 <5mass)- Thus one 
can compute the various biasing parameters in an arbitrary cosmological model at a given redshift 
and a smoothing length in a straightforward fashion. We derive the joint probability function 
assuming the primordial random-Gaussianity of the dark matter density field, and thus the results 
are sufficiently general. 

We note here, however, that our primary purpose is to present a general formulation to predict 
biasing properties of halos, and not to make detailed predictions at this point. In fact we adopt a 
few approximations with limited validity in presenting specific examples, but this is not essential 
in our paper and the resulting predictions can be improved in a straightforward manner if other 
analytical/numerical approximations become available. Nevertheless we would like to emphasize 
that our simple analytical prescription largely explains the basic features of the biasing parameters 
reported in the previous numerical simulations (Kravtsov & Klypin 1999; Somerville et al. 2000). 
Thus our prescription is supposed to capture the most important processes in the halo biasing. 

We organize the paper as follows. In §2, we describe a general formalism for the one-point 
statistics of the galaxy and the mass distributions from a point of view of the hidden variable 
interpretation of the nonlinear and stochastic biasing theory. Applying this general formalism, 
we develop a model for dark matter halo biasing treating the halo mass and formation epoch 
as the hidden variables in §3. The resulting expression for the conditional probability function, 
-P('^haiol'^mass); Can be numerically evaluated using the extended Press-Schechter theory, and we 
show various predictions for the halo biasing in §4. In particular, we pay attention to their scale- 
dependence and redshift evolution, and compare our model predictions to previous simulation 
results. Finally section 5 is devoted to the conclusions and discussion. 
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2. Formulation of Nonlinear Stochastic Biasing in Terms of Hidden Variables 

In this section, we present a general formalism of biasing for the one-point statistics of the 
galaxies and the mass smoothed with the radius R. While we specifically focus on the second order 
statistics and discuss their nonlinearity and stochasticity in the present paper, the formalism is 
readily applicable to higher-order statistics. 



2.1. Probability distribution function and hidden variables 

Recall that the fluctuations of galaxies and the dark matter density field are given by 



(^gallX, Z) = 1, (^mass(x, z) = = 1, (4) 

'^gal Pmass 

where the variables with over-bar denote the homogeneous mean over the entire universe. We 
evaluate these quantities smoothed with a spherical symmetric filter function T^(x; R): 

5(x, z\R) = J £y W{\^ - y\;R)6{y, z), (5) 

which corresponds to quantities in equation (1). The one-point statistics of galaxies and the mass 
are calculated from equation (5). 

In general, fluctuations of the biased objects Jgai are specified by multi-variate functions of 
(5inass and other observable and unobservable variables, U, O, characterizing the sample of objects. 
Then one can formally write 

,5gal = Ag(i?,z|,5^ass, U, O), (6) 

where we use Ag to denote the galaxy number density contrast at a position x and redshift z 

smoothed over a scale i? as a function of 5mass, and O. In the above expression, O and V( should 
be also regarded as functions of (x, z) smoothed over the size R. In practice, galaxies in redshift 
surveys are identified and/or classified according to their magnitude, spectral and morphological 
type. The spatial clustering of galaxies should naturally depend on those observable quantities, O. 
Since any sample of galaxies is selected over a range of values for O, the distribution of O leads to 
the stochasticity of the clustering biasing of the sample. Furthermore, the unobservable quantities 
or the hidden variables, U, which characterize an individual galaxy reflecting the different history 
of gravitational clustering and merging, radiative cooling, and environment effects, should provide 
additional stochasticity. Although these processes could be related to the dark matter density 
fluctuation in a "non-local" fashion, we intend to incorporate those effects into our biasing model 
by a set of local functions such as the gas temperature, mass of the hosting halos, and the formation 
epoch of galaxies. 

While the distinction between O and U is conceptually important, it may not be easy or 
straightforward in reality. Nevertheless it is not essential in our prescription below as long as their 
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probability distribution function (PDF): 

P = P{6m^s,U) (7) 

is specified. It sliould be noted that the above expression implicitly depends on the smoothing 
radius and the redshift for the given classes O. The statistical information of galaxy biasing is 
obtained by averaging over the joint PDF P((5mass, <^gai)- To be more specific, the joint average of 
a function J^{5jaass^ f^gai) is defined by 

mass) 

— j ^^^^^ j 

where we use (••■)) so as to expUcitly denote the joint average over the two stochastic variables, 

<^mass S'lld <5ga,l. 

In our prescription, however, it is more convenient to perform the averaging over P((5inass) U) 
instead of P((5i„ass, (^gai): 

((■^((^mass, (^gal))) = 

i 

where the variable Ag in the argument of should be regarded as a function of (^mass and U (eq.[6]). 
Of course the two expressions (8) and (9) should give the identical result, and thus one obtains 

PiSuiass,Sga\)dSga,i= ■■■ Y[ ^ i^^^s , l^) , (10) 

J JC{U) %^ 

where the region C{U) of integration is defined as 

C{U) = {U\ (5gal < Ag < ,5gal + d(5gal}, (11) 

for a given (^mass- Equation (10) implies that the unobservable information represented by U serves 
as a source for stochasticity between Smass and (5gai . In other words, equations (10) and (11) can be 
regarded as to the definition of the joint PDF P(5mass) <^gai)- Once P{SinassM) is specified, equation 
(10) can be computed numerically in a straightforward manner. In the next section, we will develop 
a simple analytical model for the dark halo biasing and explicitly calculate the joint PDF according 
to this prescription assuming that the formation epoch of halos Zf is the major variable in U. 

Finally when the joint PDF is given, it is straightforward to calculate the conditional PDF of 
galaxies for a given over density (^mass- 

P(<5gal|<5^ass) = ''^p7r'V^ ' (12) 

where the one-point PDF of the mass (^mass is related to 

(^^gal-P(<^mass,<^gal)- (13) 
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2.2. Second-order statistics and the biasing parameters 

Turn next to the statistical quantities characterizing the PDF. For this purpose, we consider 
the second-order statistics, variance of galaxies, variance of mass, and covariance of galaxies and 
mass, which are defined ^ respectively as 

agg(i?,z) = ((5gal^)), (T^^(i?,z) = ((5mass^)), Cr|m ^) = ('^gal'^mass) • (14) 

It should be also noted here that the last quantity, covariance of galaxies and mass, is not positive 
definite unlike the first two. As we will see below, however, this is always positive for a range of 

1 /2 

model parameters of cosmological interest which we survey. Thus we use (Tgm(-R, z) = ((5gai<^mass)) 
for the notational convenience. Their ratios represent the degree of the biasing and stochasticity: 



2 

"var — ) ' corr — > V-'-'-'/ 



which are sometimes quoted as "the" biasing parameter and the cross correlation coefficient (Dekel 
h Lahav 1999). In this paper, we use the subscripts, var and corr, for the above parameters in 
order to avoid possible confusions with other parameters introduced in previous papers. 

The above definitions of 6var and rcorr do not yet fully distinguish the nonlinear and stochastic 
nature of the biasing in a clear manner. Thus we introduce more convenient statistical measures, 
e„i and egcatt) which quantify the two effects separately. For this purpose, the conditional mean of 
5gai for a given (5inass (Dekel & Lahav 1999): 

<5gal(<5mass) = ^ dSga,\ PiSga.\\Smass) ^gal (16) 

plays a key role. Note that the average of 5gai('^mass) over 5mass vanishes from definition (4) : 



(^f^mass -f ('^mass) '^gal('^mass) — J" (^^gal 



d^mass ('^mass) '^gal) 



eal 



= J ddg^iPidg^i)dg^, = 0. (17) 
The nonlincarity of biasing refers to the departure from the linear proportional relation between 



5gai and (5mass- This can be best quantified by the following measure: 



^^2 ^ jOmas/} (((^galO) _ ^ ^ ^mm ('^gal )) _ ^ ^^g^ 
((<^mass<5gal))^ "^gm 



The second equality in the above comes from the fact that ((^mass'^gal)) = f'^mass'^gal))- Prom the 
Schwartz inequality, one show that the right-hand-side of the above equation is non-negative and 
vanishes only if the linear coefficient bi = (^gai/(^mass is independent of 5mass- 



^While we use ( • • ■ ) (joint average over (Jmass and Jgai) even in the definition of crgg and cTmm, it can be replaced 
by the single average over Sg^i and Jmass, respectively. 
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The stochasticity of biasing corresponds to the scatter or dispersion of 5gai around its condi- 
tional mean (^gai(<^mass)- Averaging this scatter over ^mass with proper normahzation, we define the 
following measure for the stochasticity of biasing: 

2 _ (('^mass')) (((^gal - M") _ ^mm ^gg " ..q. 

\\Omass"gal// gm 

Since the galaxy density field (5gai (eq.[6]) depends on many variables other than ^mass, ^scatt does 
not vanish in general. In turn, Cgcatt = corresponds to the unlikely case that Sg^i is uniquely 
determined by (^mass and thus 5gai = (5gai(5mass)- 

The galaxy biasing still exists even when Cgcatt = ^ni = 0. In fact, a simple linear and deter- 
ministic biasing (2) falls into this category. This effect can be separated out from the covariance 
or linear regression of Jgai and ^mass (Dekel & Lahav 1999) as follows: 

, _ f'^mass'^gal)) _ '^gm tOCi^ 

~ (5 ^)) ~ ■ ^ ^ 

W^mass // '^mm 

This quantity is equivalent to the coefficient of the leading order in the Taylor expansion, 6ggj = 
fecov<^mass + ' ' 'j in a perturbative regime, S^^ss 

< 1 (Taruya & Soda 1999). 

The biasing parameters that we introduced are related to the more conventional biasing coef- 
ficients (eq.[15]) as 

^var = ^cov(l "I" ^scatt "I" ^nl ) ^ i ^corr — / (^1) 

V 1 + ^scatt'' + Cnr 

These relations clearly indicate that egcatt and eni separate the stochastic and nonlinear effects which 
are somewhat degenerate in the definitions of ?)var and rcorr- Also it may be useful to express our 
biasing parameters in terms of those introduced by Dekel & Lahav (1999) : 



0, Cscatt — £nl — A 




(22) 



Of course the two sets of choice are essentially equivalent, we hope that our notation characterizes 
the physical meaning of nonlinear and stochastic biasing more clearly. Finally, while the present 
paper is focused on the analyses of the second-order statistics, it is fairly straightforward to extend 
the above formalism to the higher-order statistics. 



3. An Analytic Model of Nonlinear Stochastic Biasing for Dark Halos 
from the Formation Epoch Distribution 

3.1. Schematic picture of our halo biasing model 

The previous section describes the general formalism for the nonlinear stochastic biasing in 
terms of the hidden variable interpretation. In this section we present a specific model for halo 
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biasing and discuss its predictions according to the general formalism. Before proceeding to the 
technical details, it is useful to explain first the basic picture of our model in a qualitative manner. 

As illustrated schematically in Figure 1, we consider the mass and galaxy density fields at 
redshift z smoothed over the top-hat Eulerian proper radius R. The mass density contrast (^mass 
computed in the Eulerian coordinate relates R with its Lagrangian coordinate counterpart Rq = 
(1 + '-'mass 

assuming the spherical collapse. Then the mass in the sphere Mq is simply given 
by (47r/3)pmass-Ro' where Pmass(^) is the physical mass density at z. Also the linearly extrapolated 
mass density contrast Sq in the sphere can be evaluated from (^magg on the basis of the nonlinear 
spherical collapse model (e.g., Mo & White 1996). 

Each sphere of the Eulerian radius R should contain a number of gravitationally virialized 
objects i.e., dark matter halos. Given Mq and Sq, their conditional mass function can be predicted by 
extended Press-Schechter theory (e.g., Bower 1991; Bond et al. 1991). Such halos are conventionally 
characterized by their mass Mi and linearly extrapolated mass density contrast 6i assuming that 
their formation epoch is equivalent to the current redshift z. Kitayama &: Suto (1996a) pointed out 
that this approximation often significantly changes the predictions for X-ray cluster abundances on 
the basis of the Press-Schechter theory, and proposed a phenomenological prescription to compute 
the formation epoch Zf. In fact, the halo biasing derived by Mo & White (1996) is fairly sensitive 
to the difference of z and Z{ as noticed by Kravtsov & Klypin (1999). Thus we extend the biasing 
model of Mo and White (1996), in which Zf should be specified a priori, by considering explicitly 
the dependence on zj and averaging according to the formation epoch distribution function of Lacey 
Sz Cole (1993) and Kitayama &; Suto (1996b). This is a major and important improvement of our 
model over the original proposal of Mo and White (1996). In our model, therefore, the formation 
epoch Zi and the mass of halo Mi constitute the hidden variables U (§2), and their PDFs generate 
the nonlinear and stochastic behavior in the resultant halo biasing. 



3.2. Halo biasing from the extended Press-Schechter theory 



In what follows, we assume that the primordial mass density field obeys the random - Gaussian 
statistics (e.g., Bardeen et al. 1986). In this case, the (unconditional) mass function of dark halos 
(Press &: Schechter 1974) : 



n(Mi,z; Si)dMi 



1 Pi] 



27r Ml a3(Mi,z)^^^[ 2a^{Mi,z)_ 



Si 



da^{Mi,z) 



dMi 



dM^ 



(23) 



proves to be in reasonable agreement with results from A'^-body simulations (e.g., Efstathiou et al. 
1988; Suto 1993; Lacey &: Cole 1993,1994). Equation (23) corresponds to the comoving number 
density of halos exceeding the critical density threshold 5i and of mass between Mi and Mi + dMi . 
The value for Si will be specified when we consider the one-point function or the conditional PDF 
of the dark halos below (see eqs. [27] and [34]). The mass variance a^{Mi,z) is defined from the 
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linear power spectrum of mass density fluctuations Piin(fc) at present {z = 0): 

a\M,,z) = D\z) J ^ k^PUk)W\kRi), (24) 

where D{z) is the hnear growth factor normahzed to unity at z = and the top-hat window 
function: 

3 

W (x) = ^ {sin X — X cos x) (25) 

X"^ 

with the Lagrangian radius Ri = (47r/0inass/3Mi)'"^/^ is adopted. 

Since our one-point statistics of halos is evaluated within a sample of spheres of the Eulerian 
radius R, we need the conditional mass function for halos within a sphere. For this purpose, we 
use the extended Press-Schechter theory which predicts the conditional mass function for halos of 
(Ml, 5i) and (Mi -|- dMi, 5i) in the background region of (Mq, Sq) as 

n{Mi,6i\Mo,6o;z)dM, = ^ ^'""^^ 



X exp 



2^ Ml [a2(Mi,z)-a2(Mo,z)p/2 

" (^1 - dof 1 da\Mi,z) 

2{a2(Mi,z)-a2(Mo,z)} 



dMi. (26) 



Then the biased density field for halos of mass Mi, which formed at Z{ and are observed at z 
within the sampling sphere of {R, Smass) in Eulerian coordinates, is derived by Mo &; White (1996): 

('halo = ^h{R, Z\ dmass, Zf, Mi) = (1 + Vass) 7TJ FT l" i^') 

n[Mi,z;dc{z,Zi)) 
In the above, the critical threshold Sc for those halos is given as 

Sc{z, z,) = 4,0 5c,o - ^^^^ - 1.69, (28) 

again on the basis of the spherical collapse model. The remaining task is to compute the mass, Mq, 
and the linearly extrapolated mass density contrast, 5q, of the sampling sphere from its Eulerian 
radius and density contrast (-R, (5mass)- Since 5q^1, the former is simply given as 

47r q 47r q 

Mo = Y Pmass(^) (1 + 5o) ^0 - y Pmass(^) (1 + (^mass) , (29) 

with pui3ss{z) being the physical mass density at z. Finally we adopt the following fitting formula 
obtained in the spherical collapse model (Mo & While 1996) to compute 5q in terms of S^ass- 

1.35 0.78785 1.12431 

° " (1 + s^^^yi^ + (1 + <5^a.s)°-5«6«^ (1 + s^^y/^ ■ ^-^^^ 

While equations (28) and (30) were originally derived in the Einstein-de Sitter model, we numeri- 
cally checked that they provide sufficiently accurate approximations for our present purpose even 
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in the rio = 0.3 and Aq = 0.7 model that we consider later. Thus we use the expressions (28) and 
(30) irrespectively of the cosmological models in the subsequent analysis. 

Figure 2 illustrates the dependence of Ah on M and Zf as a function of (Jmass smoothed over 
R = 8/i~^Mpc at z = 0. Given a halo mass, Ah is very sensitive to the formation epoch Z{ especially 
in the range of Zf ^ z + 1. As Zf increases, the dependence Ah and M and (Jmass becomes significant; 
this reflects the fact that the larger mass halos preferentially form in the denser regions than the 
average since the typical halo mass that can be collapsed and virialized decreases at higher . Such 
M and zi dependence of Ah convolved with the PDF of M and Zf leads to the nonlinear stochastic 
behavior of the biasing of dark halos. 

While we regard halo mass Mi and its formation epoch Zi as the hidden variables in equation 
(27) , they may not be entirely unobservable. One may infer the halo mass and the formation epoch 

for an individual galaxy by combing the observed luminosity, color and mctallicity with a galaxy 
evolution model. In this case, their probability distribution functions need to be convolved with 
such observational selection functions with our prior distribution. Except for this correction, our 
methodology presented below remains the same. 



3.3. Formation epoch distribution of dark halos 

As indicated in equation (27) , the amplitude of halo biasing is explicitly dependent on its for- 
mation epoch Z{. Thus the simple approximation of Zf = z may lead to even qualitatively incorrect 
predictions for the biasing. In fact this was shown to be the case in recent N-body simulations 
by Kravtsov & Klypin (1999). The importance of the distribution of Zf has been emphasized by 
Kitayama & Suto (1996a) in a different context, and a model for its PDF was proposed by Lacey 
&; Cole (1993). Incidentally Catelan et al. (1998a,b) also proposed a different model of halo bi- 
asing considering the Zf-dependence. Their model simply treats Zf as a free parameter and does 
not properly take account of its distribution function. Our model presented here incorporates the 
distribution function of the formation epoch explicitly. 

Adopting the excursion set approach (Bond et al. 1991) and defining the formation redshift 
Zf of a particular halo of mass M at z as the epoch when the mass of its most massive progenitor 
exceeds M/2 for the first time, Lacey Sz Cole (1993) derived the differential distribution of the halo 
formation epoch dp/dz{. Their result is expressed as 

g(.HM,.)cf.f = ^(^f,M)g.., (31) 

M) = ' r 4 f 1 - 3) e-.V(2^), (32) 

dui^ ' V2^Jo 53/2M'(5)V Sj ^ ' 

where 

Sc{z, Z{) - (5c,o 



ibi{z{,M,z) = 



^ya'^{M/2,z)-a^{M,zy 
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S{M',M) = 



a 



a^{M',z) -a^M,z) 
72(M/2,z) -(72(M, z)' 



(33) 



The function M'{S) in the integrand of equation (32) can be obtained from equation (33). While 



the above expressions are rather compHcated, practical fitting formulae for the mass variance in cold 
dark matter (CDM) models and for the formation epoch distribution were obtained in Kitayama 
&; Suto (1996b) which we adopt throughout the analysis. Those are summarized in Appendices A 
and B for convenience. 

It should be noted that the definition of halo formation is somewhat ambiguous in the frame- 
work of the extended PS theory. This aspect is explored in Kitayama & Suto (1996a), and their 
Figure 1 explicitly shows how the result is dependent on the adopted ratio of the current halo mass 
and the progenitor mass at the formation epoch. The figure implies that the resulting formation 
rate is fairly insensitive to the value around 0.5 that we adopt here. 

Figure 3 plots the formation epoch distribution for halos selected at z in CDM models. 

Specifically wc choose {Qo,Xo,as,h) = (0.7,0.3,1.0,0.7) (Lambda CDM; hereafter LCDM ) and 
(1.0,0.0,0.6,0.5) (Standard CDM; hereafter SCDM). The top-hat mass fluctuation amplitude at 
8h~^yipc, as, is normalized to the cluster abundance (Kitayama & Suto 1997). We show results 
for halos of M = lO^^/i-^M© (dashed), 10^^Mq{soM) and lO^^/i-iM© (dot-dashed) in LCDM, and 
10^2j\fQ((jotted) in SCDM. The shape of those distribution functions is quite similar, and charac- 
terized by a peak around Zf = z + (0.5 ~ 1). The peak redshift becomes closer to the observed 
one, z, and the distribution around the peak becomes narrower as the halo mass increases, both of 
which are easily understood in the hierarchical clustering picture like the present models. 

Note that the SCDM model generally predicts a more sharply peaked distribution closer to z 
than the LCDM model (compare solid and dotted lines in Fig. 3). This is also reasonable from the 
fact that the growth of fluctuations is rapid in SCDM and thus halos form only recently. Thus the 
formation epoch distribution is fairly sensitive to the cosmological parameters. 

3.4. Conditional and joint probability distribution functions of dark halos and mass 

Now we are in a position to explicitly construct the conditional probability distribution of 
the dark halo -P((5haio|<^mass) for a given (Jmass; and the joint probability distribution P{Sha.io, <5mass)- 
Basically we follow the prescription described in §2.1, but in slightly different order. 

We first compute -P(5haio|'^mass) applying equation (10): 




(34) 



where the region C(M,Zf) of the integration is determined from the following conditions: 



C{M,Zi) = { {M,Zi) I <5halo < Ah(i?,z|<5, 



M,Zi) < 4alo + C?<5halo, 
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Minin <M< M^^, Z<Zi<(x} (35) 

The normalization factor J\f is defined as 

M= dM dzi -^{zi\M,z) n{M,z;6c,o). (36) 

The integrand in equation (34) just corresponds to the joint PDF P{U\dinass) for a given 5mass- 
Since the joint PDF is simply computed according to 

-f ('^mass) '^halo) = -P('^halo|<^mass) -P('^mass)) (37) 

one needs a reliable model for the one-point PDF of dark matter density contrast, P(5mass)- Fortu- 
nately it has been known that this can be empirically approximated by the log-normal distribution 
function to a good accuracy (e.g., Coles & Jones 1991; Kofman et al. 1994; Bernardeau & Kofman 
1995; Taylor & Watts 2000): 



-f ('^mass) ^'^mass — n- — exp 



2al 



1+6, 



where 

al = Hl + alj, (39) 

and cJmm is defined in equation (14). Note that equation (38) reduces to the Gaussian distribution 
for (Tmm ^ 1, and thus this model again assumes the primordial random-Gaussian density field 
implicitly as our entire analysis. Finally we adopt the fitting formula (Peacock & Dodds 1996) for 
the nonlinear CDM power spectrum P^iik, z) in computing the mass variance: 



2 



{R,z) = / H k^Pni{k,z)W\kR) , (40) 



with W{kR) being the top-hat smoothing function (eq.[25]). The validity of the log-normal approx- 
imation for the one-point PDF is examined by Bernardeau Sz Kofman (1995); their Figure 10 indi- 
cates that equation (38) reproduces the simulation results very accurately at least for R ^ 5/i~^Mpc. 
Although the accuracy on smaller scales is not shown quantitatively, it would be reasonable to as- 
sume that the approximation is acceptable up to i? '-^ l/i~^Mpc. Also our statistical results are not 
sensitive to the tail of such PDF in any case. 

Substituting the analytical expressions for dp/dz{{zf\M, z) and n{M, z; 5cfi) discussed in §3.3 
into equation (34), one may numerically compute the conditional PDF i-*((5halo|'5mass)) and the joint 
PDF P(5haio, <^mass) from equation (37). Then all the statistical quantities can be evaluated using 
equation (8). In practice, however, it is more convenient and even accurate to use (9) which in the 
present case is written explicitly as 

CO 



/CO 
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/ dM dzf ^(Ah, 

J Mm\r. J Z 



6„ 



dp 
dzf 



{zf\M,z) n{M,z;6c,o) 



■ (41) 



We use the above expression in evaluating the various biasing parameters below except in presenting 
the PDFs directly. 



4. Results in CDM models 

Wc present several specific predictions applying our nonlinear stochastic halo biasing model to 
representative CDM models mentioned in §3.3. Throughout the analyses, we consider the range of 
halo mass between Mmin = lO^^hr^MQ and M^ax = W^^h'^MQ unless other wise stated. 

The general formulation described in the previous section should work in principle even on fully 
nonlinear scales. In practice, however, the results presented below are limited by the halo exclusion 
effect (due to the finite size of halos) and our approximation, equation (38), for the one-point PDF. 
The validity of both effects should be carefully checked on small scales. Since a typical virial radius 
of a halo of mass M in LCDM model is ~ 0.5{M/10^^ MQ)^/^h~^ Mpc, the exclusion effect cannot 
be neglected below R = l/i~^Mpc but is not so strong for R ^ 3/i~^Mpc even for our largest mass 
considered (Mmax = Mq). Also the validity of of the log-normal approximation is already 

remarked in subsection 3.4. Thus we expect that our predictions below are fairly reliable up to 
R ~ 3/i-^Mpc. 



4.1. Conditional probability distribution of ^mass and S^aXo 

Since the conditional PDF P(5haio|'^mass) plays a central role in the Dekel & Lahav (1999) de- 
scription of the nonlinear stochastic biasing, we first present -P((5haio|<^mass) predicted from our 
model. For this purpose, we start with equations (34) and (35). Specifically we divide the 
(•^gah <^mass) plauc in a 3000 X 3000 mesh, and accumulate the integrand of equation (34) satis- 
fying the constraint (35) on each grid. The resulting PDFs are plotted in Figure 4 for a given 
mass density contrast; Jmass = (solid), d^aass = 1 (dotted), 5mass = 2 (dashed) and (Jmass = 3 
(dot-dashed). The upper and lowers panels show the results at ^; = and z = 1, respectively, with 
the top-hat smoothing radius of i? = 8^~^Mpc {left) and R = 20^~^Mpc (right). 

The ticks on the upper axis in each panel indicate the corresponding conditional mean (^haio 
(eq.[16]). The peak position of the distribution is in reasonable agreement with (5haio- As Figure 
2 indicates. Ah given z and R is fairly monotonically dependent on M and Zf. Thus the peak in 
the conditional PDF reflects that of the formation epoch distribution dp/dz{. The width of the 
distribution around the peak, on the other hand, is dominated by the mass distribution since Ah 
becomes more sensitive to the halo mass in a denser environment (Fig. 2). 
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4.2. Joint probability distribution of (Jmass and S^aio 

Once -P((^haio|^mass) is given, the joint PDF P(5mass, (^haio) is simply obtained by multiplying 
the one-point PDF of the mass density P{6mass), in our case, the log-normal model (eq.[38]). The 
resulting contour on (^mass, <^haio) plane is illustrated in Figure 5. This example shows the result 
with the top-hat smoothing radius of i? = 8/i""^Mpc at different redshifts. Solid lines in each panel 
indicate the conditional mean 5haio(<^mass)- 

The number density of halos of mass exceeding the current threshold Mmin = lO^^/t^^M© 
become progressively smaller as z increases. Such halos naturally reside in higher density regions, 
and therefore are strongly biased with respect to mass. The biasing of those halos gradually 
decreases as time since they simply follow the gravitational field of the background mass after 
formation (Fry 1996; Tegmark &; Peebles 1998; Taruya, Koyama & Soda 1999). In addition, 
new halos with M > Mmin form more easily later and can be found even at moderately dense 
environment. For both reasons, the mean bias t^haio as a function of S^ass decreases at lower 
redshifts. 

At z = 0, the joint PDF -P(5mass) ^haio) shows slightly anti-biasing behavior, i.e, (5haio ~ <^mass- 
This is partly due to the fact that a fraction of halos with M < Mmax are merged into a part of 
larger mass halos since the typical virialized halo mass M*(z) at z = approaches the mass scale 
Mo{R), corresponding to our adopted smoothing radius R = 8/i~^Mpc itself. In other words, our 
halo model generally predicts the positive-biasing except for those halos of M ^ M*(z) ^ Mq{R) 
for given R and z. 

While Figure 5 elucidates the global feature of the joint PDF, the statistical weights are 
practically dominated by the relatively narrow regions around 5haio ~ <^haio(<^mass)- Those regions are 
illustrated better in logarithmic scales rather than in linear scales. For this purpose, we recompute 

the joint PDF from equations (34) and (35) using the 3000 x 3000 mesh with the logarithmically 
equal bin on the (1 + (5gai, 1 + 5mass) plane. Note that the resulting PDF sampled in this way, 
-P(l + ^mass, 1 + <^haio), Satisfies 

-P(l + (^inass, 1 + (^halo) C?ln(l + (5mass) 6^ ln(l + (^halo) = ^'('^mass, ^halo) C^'^mass C^^halo- (42) 

Thus we decide to plot the joint PDF (l+(5mass)(l+'^haio) -P(<^mass; <5haio)- In Figure 6, the dispersion 
around the mean biasing decreases at higher z and/or larger R in contrast to the conditional PDF 
plotted in Figure 4. This is basically because the PDF P((5niass) in linear regime becomes toward 
Gaussian and sharply peaked. 

We remark that the contours of our joint PDF plotted in Figures 5 and 6 seem to be very 
narrow around ^mass ~ (^halo ~ 0. This can be understood by the fact that the positive halo density 
contrast preferentially developed on the over-dense dark matter environment. In other words, this 
is a natural consequence of our bias model in which the signs of (5haio and 5mass are almost the same 
as illustrated in Figure 2. Incidentally this feature may be visually exaggerated by the contours of 
very small probabilities; if one focuses only on the contours of P((5mass5 <^haio) > 0.01, the effect does 
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not look so strong. In reality, and thus in numerical simulations, additional stochastic processes 
other than the mass and formation epoch distribution (including the dynamical motion of halos) 
should further increase the scatter which makes the contour rounder than our predictions. 

We will return to these contour plots in understanding the behavior of the biasing parameters 
later. 



4.3. Scale-dependence 

The previous two subsections show that the joint and the conditional PDFs are dependent on 
both R and z. This dependence is translated to the scale-dependence and redshift evolution of the 
second order statistics defined in §2.2. 

Figure 7 shows dhh) CThm and cJmm at different rcdshifts as a function R. We compute those 
statistics directly integrating equation (41) over M, Zf and ^mass) instead of using -P(5mass) <^haio)- 
While their scale- and time-dependence is noticeable even from those panels, the biasing parameters 
(^var, ''corr, ^covj ^scatt ^nd plotted in Figures 8 and 9 are more suitable in understanding the 
origin of the behavior. 

Consider first the scale-dependence (Fig.8). While 6var and 6cov are generally a decreasing 
function of i?, this behavior is significant only up to i? 2(1 + 2;)/i~^Mpc in this model. This 
feature is more quantitatively exhibited by e^i and egcatt {Lower-right panel). Therefore in practice 
the linear biasing provides a good approximation on linear and quasi-linear regimes. The biasing 
is non-deterministic especially on smaller scales almost independently of z (see also Fig. 9). In 
addition, rcon does not approach unity even on large scales, implying that non-deterministic nature 
still exists there to some extent. As the lower-right panel in Figure 8 indicates, egcatt S> eni on all 
scales and the above feature should be ascribed to the stochasticity due to the distribution of M 
and Zf . In fact, this stochasticity on large scales is expressed explicitly in terms of the linear biasing 
approximation by Mo &: White (1996): 

z/2 - 1 

5halo - ^'mw <^mass, hyi^{z\M , Zf) = 1 + — -, (43) 

dc[Z, Zi) 

with v = Sc{z, Zi)/a{M, z). It should be noted that 6mw is often regarded as a function of z and M 
assuming z = Zf, leading to the linear deterministic model. Thus once a halo mass M is specified, 
Cscatt = 0. In equation (43), however, we explicitly keep the Zf-dependcncc which adds the stochastic 
nature in the model. More specifically, the definition (19) with equations (41) and (43) reduces to 

^ 2 _ (('^mass^)) ((^MW<^mass " ((&Mw))M,2f '^mass)^)) _ i^Mw) M,Zf ^^^^ 
(((^Mw))M,2f<^mass^)) (^Mw))M,Zf 

where {■ ■ ■}m zt denotes the average over M and Z{. This accounts for the scale-independent non- 
vanishing Cscatt exhibits in Figure 8. In conclusion, our model implies that the halo biasing does 
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not become fully deterministic even on large scales where its nonlinearity is negligible. This result 
is not surprising since we take into account the stochastic processes which do not vanish on large 
scales, but the overall effect is quite small ( 0.98). 

4.4. Redshift dependence 

Next discuss the redshift dependence of our biasing model. Figure 9 shows that 6var and 
6cov strongly evolve in time. In fact, this is in marked contrast with predictions on the basis 
of phenomenological linear deterministic models. For example, a model of Fry (1996) leads to 
evolution of a form: 

b{z) = l + ^^[b{z = 0)-l]. (45) 

This implicitly assumes that all the objects of interest form at the same Zf and that their biasing 
parameters at Z{ are independent of the mass. Since neither the above assumptions apply to our 
model, the prediction (45) is quite different from ours even in the linear regime. Our results 
generally show much stronger evolution as z increases despite the fact that b{z = 0) is very close to 
unity. The recent compilation of the various galaxy catalogs also indicates that the prediction (45) 
does not reasonably describe the behavior at z ^ 2 (Magliocchctti ct al. 1999). Thus, the proper 
modeling in the framework of the nonlinear stochastic biasing is important even in predicting 6var 
and 6cov 

In our halo biasing model, the degree of stochasticity egcatt is almost constant in time because 
it is determined by the effective widths of the probability distribution functions of M and Z{. 
Incidentally the nonlinearity eni does not evolve monotonically (thin lines in the Lower-right panel). 
At an intermediate redshift, Cni reaches at a minimum. This behavior is qualitatively explained 
from the curvature of the conditional mean (5haio as a function of (5mass, i-e., its second derivative 

mass 

\ At 

(^mass ^ 1) halos of the mass -/Vfj-Qin < M < -Mmax that we adopt here exhibit 
stronger positive biasing ((5haio > ^mass) on average at z ^ 2 and mildly anti-biasing ((5haio ~ (^mass) 
at z ^ 1. Since 5haio = — 1 at t^mass = — 1 by definition (cf., eq.[27]), the dependence results in 
positive and negative curvatures of 5haio('^mass)5 respectively at 2; ^ 2 and ^ ^ 1, especially around 
<^mass ~ where the contribution of the joint PDF -P((5mass) <^haio) is significant. This feature is 
clearly visible in Figure 5. In turn, the curvature should be minimum somewhere between z = I 
and 2. When higher-order correction is neglected, eni is dominated by the curvature or the second 
derivative of 5haio(<^mass) properly averaged over Smass (eq.[18]), and it should become minimum at 
the same redshift. It is interesting to note that the qualitatively similar evolutionary feature was 
found from the numerical simulations of Somerville et al. (2000; their Fig.17). 
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4.5. Origin of stochasticity and cosmological model dependence 

So far we have presented model predictions for halos averaged over lQ^^h~^MQ < M < 
M^^H^^Mq in LCDM model taking into account the appropriate Z{ distribution. Figures 10 and 
11 compare those fiducial results with predictions based on the different model assumptions. 

First we address the question of the origin of the stochasticity. Since our biasing model has two 
hidden parameters, M and Z{, we attempt to separate the two sources by fixing M = lO^^h^^MQ or 
Zf = z + 1 while keeping the other parameters exactly the same. The upper-panels and the lower- 
left panel of Figure 10 suggest that the Z{ distribution dominates the stochasticity at low redshift 
(z = 0), while the effect of the M distribution becomes significant at higher redshift (z = 3). The 
same behaviors can be seen in the lower-right panel of Figure 11. Since our model relies on the 
hierarchical picture of structure formation, the result is simply deduced from the merging history of 
halos. Thus, in general, the major contribution to the joint PDF can become the formation epoch 
distribution. This is also indicated in the scale-dependence of the stochasticity in the lower-left 
panel of Figure 11 (thick-dashed and thick-dotted lines). 

Next consider the cosmological model dependence. For this purpose, we plot the result in the 
SCDM model with the same mass range. The joint PDF at z = and R = 8/7"iMpc (Lower- 
right panel of Fig. 10) is confined in a slightly narrower regime compared with that in LCDM. 
This comes from the fact that the formation epoch in SCDM shows a bit more sharply peaked 
distribution dp/dz{ in than that in LCDM with the same halo mass M (Fig.3). As a result, the 
stochasticity in SCDM is smaller (i.e., egcatt is smaller and rcorr is closer to unity), but 6var at z = 
is almost insensitive to such small changes. Rather the major difference between LCDM and SCDM 
is the redshift evolution of &var which increases more rapidly in SCDM reflecting the faster growth 
rate of density fluctuations. 

4.6. Comparison with previous results 

Dark matter halos arc quite natural and likely sites for galaxy and cluster formation. Thus 
there are many previous papers to discuss different aspects of the halo biasing on the basis of 
different assumptions and modeling (Catelan et al. 1999a,b; Blanton et al. 1999). Among others, 
Kravtsov & Klypin (1999) and Somerville et al. (2000) analyzed the nonlinearity and stochasticity 
in halo and galaxy biasing using numerical simulations. In this sense their work is complementary 
to our analytical modeling, and deserves quantitative comparison with our results. 

Kravtsov &; Klypin (1999) performed high-resolution N-body simulations employing N = 256^ 
particles in a periodic 60/i~^Mpc box so as to overcome the halo over-merging. In particular, 
their Figure 4 plotting the joint PDF i-*((5niass) ^haio) is quite relevant for the comparison with our 
Figure 5. Strictly speaking, their simulated halo catalogue is based on slightly different identifi- 
cation scheme (Klypin et al. 1999); the bound density maxima algorithm, a selected mass range 
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[-^min, -^max] IS limited both by the numerical resolution (individual particle mass) and the simula- 
tion boxsize. Furthermore their statistics is based on the smoothing length R = 5/i~^Mpc. Despite 
such quantitative difference, both contours look surprisingly similar. They pointed out that the 
Mo &: White biasing with Z{ = z -\- 1 phenomenologically fits the conditional mean 5haio(^mass) of 
their simulations. Actually this is exactly what we find here, and mainly due to the fairly strong 
peak in the formation epoch distribution (Fig. 3). 

Somerville et al. (2000) examined the nonlinear and stochastic nature in the galaxy biasing 
combining N-body simulations and semi-analytic model of galaxy formation (see also Kauffmann, 
Nusser & Steinmetz 1997). While their model incorporates many detailed astrophysical processes 
(gas cooling, star formation, stellar evolution and population synthesis, galaxy feedback and merg- 
ing among others) that our analytical modeling neglects, their overall conclusion is that the massive 
halos (M ^ IO^^/i^^Mq) identified from simulations have nearly one-to-one correspondence with 
luminous galaxies. Their conclusion is encouraging to us since it justifies our crucial assumption 
that the nonlinearity and stochasticity in the halo biasing are originated mainly from gravitational 
clustering. In fact, their Figs. 5 and 6 are also very similar to our joint PDF plotted in Figure 
6. Furthermore they find (their Fig. 17) that both nonlinearity and stochasticity of galaxies evolve 
moderately as redshift, and decrease on larger scales, and that the stochasticity shows a minimum 
at an intermediate redshift. These findings are fully consistent with our results and in fact can be 
explained physically in the framework of our analytic description. 



5. Discussion and Conclusions 

In this paper, we presented a general formalism to describe the nonlinear stochastic biasing 
from the hidden variable interpretation. According to this formulation we proposed a physical 
model for the halo biasing assuming that mass and formation epoch of dark halos act as the 
major hidden variables. In particular, this model for the first time provides an analytic expression 
for the joint probability distribution function of the halos and the dark matter density fields. 
The specific functional forms for the PDFs can be derived, or are at least consistent with the 
phenomenological results, from the standard gravitational instability theory and the assumption of 
the random-Gaussianity of the primordial density field. Therefore we expect that the basic features 
of the nonlinear and stochastic biasing predicted from our model are fairly generic. In fact, detailed 
comparison with the previous numerical simulations by Kravtsov & Klypin (1999) and Somerville et 
al. (2000) shows good agreement with our predictions, indicating that the nonlinear and stochastic 
nature of the halo biasing is essentially understood by taking account of the distribution of the halo 
mass and the formation epoch. 

Then we introduced a set of biasing parameters from second-order statistical quantities fol- 
lowing Dekel &; Lahav (1999), which properly quantify the one-point statistics of nonlinear and 
stochastic biasing. As specific examples, we explicit compute those parameters in LCDM model, 
and examined their scale- and time-dependence. Our major findings are summarized as follows: 
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1. The biasing of the variance, 6vax = V''^hh/f"mm! evolves strongly as redshift. While its scale- 
dependence becomes noticeable on small scales and/or at high redshifts, a simple linear biasing 
model provides a reasonable approximation roughly at R ^ 2(1 + z)h~^Mpc. 

2. The stochasticity, rcon- = <^hm ^ / V ^hh Cmm exhibits moderate scale-dependence especially on 
R ^ 20ft.^^Mpc, but is almost independent of z. Since egcatt ^ Cnl in general, the stochasticity 
of halo biasing is mainly generated by the scatter which is dominated by the formation 
epoch distribution at lower redshifts and by the halo mass distribution at higher redshifts. 
Importantly, the stochasticity still remains constant even on large scales, which yields the 
small deviation of rcorr from unity. 

The fact that biasing approaches linear and scale-independent on large scales is just expected from 
previous theoretical work (e.g.. Coles 1993; Scherrer &: Weinberg 1998) and was recently suggested 
observationally (Seaborne et al. 1999). Our model further implies that the scale-independence still 
holds even in the quasi-linear regime, and more importantly, that biasing is still stochastic to some 
extent on those scales. 

Before closing we would like to comment on several limitations of our model and on future 
prospects in order. First, our analysis is entirely dependent on the formation redshift distribution 
derived by Lacey & Cole (1993). It is known, however, that the distribution becomes negative in 
some models, suggesting that this function is not fully correct. Another proposal for the distribution 
function by Kitayama & Suto (1996a) also suffers from a different conceptual problem. This might 
be ascribed to the difficulty of distinguishing the merger and accretion unambiguously within the 
framework of the extended PS theory. Despite this limitation, the two proposals by Lacey &; Cole 
(1993) and by Kitayama & Suto (1996a) lead to very similar predictions and thus we believe it 
unlikely that this significantly affects the final results. Second, several authors have claimed that a 
better agreement with numerical simulations is attained by incorporating the non-spherical effect 
into the Press-Schechter theory (Lee & Shandarin 1998; Sheth, Mo, &: Tormen 1999). Although 
we do not think that the non-spherical effect drastically changes the statistical predictions of our 
model, this is definitely a well-defined and interesting possibility to improve our model provided 
that the corresponding distribution function of the formation epoch can be derived. Third, for 
the proper comparison with observation, our predictions should be corrected for the redshift-space 
distortion. Since the distortion induced by peculiar velocity field also leads to the stochasticity, the 
resulting biasing becomes object-dependent (Taruya et al. 2000). Finally and more importantly, 
the astrophysical effects other than gravity must be taken into account. Since a halo identified in 
our scheme carries its formation epoch, it is not so diflficult to combine with the phenomenological 
galaxy formation model. Such an improved model will enable more generic and testable predictions 
for the galaxy biasing. 
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Appendices 



A. Mass variance for the CDM model 

The mass variance a{M, z = 0) defined by equation (24) requires the hnear power spectrum 
of density fluctuations Piin{k). The fitting form of the CDM power spectrum is given by Bardeen 
et al. (1986) with the scale-invariant Harrison-Zel'dovich initial condition: 

2 

[1 + 3.89g + (16.1g)2 + {5A6qf + (e.Tl?)^] "^^^ , (Al) 

where q = k/{Th Mpc~^). Here, we adopt F = ^lQhexp{—Q\) — ^/2h^l\)/QQ) with the baryon density 
parameter f^b = 0.02h~^. 

Numerical integration of equation (24) is straightforward, but time-consuming since we heavily 
use as well as its derivative da'^/dM. Fortunately Kitayama &: Suto (1996b) obtained the 
following accurate fitting formula for whose derivative simultaneously fits da'^/dM: 

(7 oc (1 + 2.208mP - 0.7668m^P + 0.7949m^P) ^^^^^^^ , (A2) 

where p = 0.0873, and m = M{Th)'^ / {10^'^ Mq). The above approximation holds within a few 
percent for both and da"^ /dM in the range 10^^ ^ m ^ 10^. 

The normalization of is characterized by the parameter as- 

a{RM = Sh'^Mpc, z = 0) = as, (A3) 

where Ru = [3M/(47rpmass)]^^^- Throughout the paper we adopt the above formula (A2) combined 
with the cluster normalization for erg (Kitayama &; Suto 1997). 



P\mik) oc k 



ln(l + 2Mq) 
2Mq 



B. Fitting formulae for the distribution function of halo formation epoch 



The distribution function of the halo formation epoch (eq.[32]) plays a central role in our 
model, but it requires a time-consuming numerical integration and inversion. Thus in the present 
paper we use the following fitting formulae of Kitayama & Suto (1996b): 

Aia) 



dp_ _ ^ 



_g-5a;f + 2C{a)ui erfc 



LUf_ 



where erfc(a:) is the complementary error function, and 



A{a) = -\/-(l-a)(0.0107-F 0.0163a), 



B{a) = 



Aia) 



C{a) 



2° - 1" 



a 



C{a) = 1 



1 — a 
25 ■ 



(Bl) 

(B2) 
(B3) 
(B4) 
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The parameter a is related to the spectral index of the mass variance (t(M). Kitayama & Suto 
(1996b) showed that in the CDM model this parameter should be replaced by 

a = Q;eff (0.6268 + 0.3058aeff), (B5) 

where the effective spectral index aeff = "C^logCcDM/'^^'^S-^ computed from the fitting formula 
(A2) and its derivative. 
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Fig. 1. — A schematic picture of our halo biasing model at a given redshift z. The mass and halo 
density fields are top-hat smoothed with the Eulerian proper radius R, and Rq is its Lagrangian 
counterpart. Each sampling sphere is characterized by R and the Eulerian mass density contrast 
Smass, or equivalently by the enclosed mass Mq and the linearly extrapolated mass density contrast 
Sq. Within each sampling sphere, there are a number of dark halos characterized by their mass 
Ml and the linearly extrapolated mass density contrast 5i. With a given z, 6i is related to the 
formation epoch Z{ of the halo. 
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Fig. 2. — The halo biasing field smoothed over R = 8/i"^Mpc at z = function of its mass 
M and formation redshift in dotted (M = lO^^/i-^M©), dashed (M = lO^^^-ij^^^)^ g-nd dot- 
dashed (M = lQ^^h~'^MQ) lines. Upper-left: Zi = 0; Upper-right: Zi = 0.5; Lower-left: Zi = 1.0; 
Lower-right: Z{ = 1.5. 
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Fig. 3. — Formation epoch distribution of dark halos. Dotted lines indicate the results for M = 
IO^^/i^^A/q in SCDM model, while dashed, solid and dot-dashed lines are respectively for M = 
lO^^h-^MQ, and W^^h^^MQ in LCDM model. Upper-left: z = 0; Upper-right: z = I; 

Lower-left: z = 2; Lower-right: z = 3. 
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Fig. 4. — Conditional probability of Jhaio for 5mass = (solid), 1 (dotted), 2 (dashed) and 3 (dot- 
dashed) in LCDM model for IO^/i-^Mq < M < IQi^^h-^MQ. Upper-left: z = and i? = 8/i-iMpc; 
Upper-right: z = and R = 20/i-iMpc; Lower-left: z = 1 and R = 8/i^^Mpc; Lower-right: z = 1 
and R = 20/i"iMpc. 
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Fig. 5. — Contour plot of the joint probability distribution of 6mass and (5haio, -P('5mass, <5haio)- These 
examples are for IO^^/j-^Mq < M < 10^^h~^MQ and R = Sh'^Mpc in LCDM model. Upper-left: 
z = 0; Upper-right: z = 1; Lower-left: z = 2; Lower-right: 2 = 3. 
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Fig. 6. — Contour plot of the joint probability distribution of f (l + (5mass; 1 + 'J halo ) = (1 + ^ mass )(1 + 
^haio)-P(<Jmass,'5haio)- Tliese examples are for IO^^/i^^Mq < M < lO^^/i^^Afo in LCDM model. The 
redshifts z are 0, 1 and 3 from top to bottom, and the smoothing lengths R are 3/i~^Mpc, 8/i~^Mpc 
and 20/i~^Mpc from left to right. 



- 30 - 




Fig. 7. — Scale-dependence of a^h (solid), iThm (dashed), and amm (dot-dashed) at different redshifts; 
Upper-left: z = 0; Upper-right: z = 1; Lower-left: z = 2; Lower-right: z = 3. 
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Fig. 8. — Scale-dependence of biasing parameters, 6var, "Tcorr, ^cov, Cscatt and 6^1, at the different 
redshifts z = (solid), z = 1 (dot-dashed), z = 2 (dashed), and z = 3 (dotted). In the Lower- 
right panel, thick and thin lines indicate the degree of stochasticity egcatt; and nonlinearity eni, 
respectively. 
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Fig. 9. — Evolution of biasing parameters, 6var; ^'corr, ^cov, ^scatt and eni, at i? = Sh^^Mpc (dotted), 
(dashed), and (solid). In the Lower-right panel, thick and thin lines indicate 

the degree of stochasticity egcatt; and nonlinearity e^i, respectively. 
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Fig. 10. — Joint probability distribution of (5mass and (^haioj -P(<5massj ^haio) at R = 8h ^Mpc for 
different models. Upper-left: LCDM with Zf = z + 1 at z = 0; Upper-right: LCDM with M = 
IO^^/i^^Mq at z = 0; Lower-left: LCDM with Zf = z + 1 and with M = lO^^ft-^^M© (the value of 
(^mass is artificially shifted by five for clarity), at z = 3; Lower-right: SCDM with W^^h"^ Mq < 
M < at z = 0. 
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Fig. 11. — Comparison of ftvan ^'corr, Escatt and eni in different models. Left: at z = 0; Right: at 
R = 8/i"^Mpc. Different lines indicate results for LCDM with M = lO^^/j-^M© (dotted), LCDM 
with Zf = z + 1 (dashed), and SCDM with lO^^h'^MQ < M < lO^^/i^^M© (dot-dashed), in 
comparison with our canonical model (LCDM with W^^ Mq < M < 10^^h~^MQ) plotted in 
solid lines. 



